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We show how the quantum to classical transition of the cosmological fluctuations produced during 
inflation can be described by means of the influence functional and the master equation. We split 
the inflaton field into the system-field (long- wavelength modes), and the environment, represented 
by its own short-wavelength modes. We compute the decoherence times for the system-field modes 
and compare them with the other time scales of the model. 
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I. INTRODUCTION 

The emergence of classical physics from quantum be- 
haviour is important for several physical phenomena in 
the early Universe. This is beyond the fundamental re- 
quirement that only after the Planck time can the metric 
of the Universe be assumed to be classical. For exam- 
ple, the inflationary era is assumed to have been induced 
by scalar inflaton fields, with simple potentials 1]. Such 
fields are typically assumed to have classical behaviour, 
although in principle a full quantum description should 
be used. The origin of large scale structure in the Uni- 
verse can be traced back to quantum fluctuations that, 
after crossing the horizon, were frozen and became clas- 
sical, stochastic, inhomogeneities 0- 

It is generally assumed that several phase transitions 
have occurred during the expansion of the Universe 
As in the case for the inflaton fields, the (scalar) or- 
der parameter fields that describe these transitions are 
described classically. However, the description of early 
universe phase transitions from first principles is intrin- 
sically quantum mechanical Q . As a specific application 
of the previous point, the very notion of topological 
defects (e.g. strings and monopoles) that characterize 
the domain structure after a finite-time transition, and 
whose presence has consequences for the early universe, 
is based on this assumption of classical behaviour for the 
order parameter 6], as it distributes itself between the 
several degenerate ground states of the ordered system. 

In previous publications, we analysed the emergence of 
a classical order parameter during a second order phase 
transition and the role of decoherence in the process of 
topological defect formation 0; H, ^| ■ 

In the present paper our concern is directly related 
with the first point above, the quantum to classical tran- 
sition of the inflaton. Any approach must take into ac- 
count both the quantum nature of the scalar field and the 
non-equilibrium aspects of the process . The problem 
of the quantum to classical transition in the context of 
inflationary models was first addressed by Guth and Pi 
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|12| . In that work, the authors used an inverted harmonic 
oscillator as a toy model to describe the early time evo- 
lution of the inflaton, starting from a Gaussian quantum 
state centered on the maximum of the potential. They 
subsequently showed that, according to Schrodinger's 
equation, the initial wave packet maintains its Gaussian 
shape (due to the linearity of the model) . Since the wave 
function is Gaussian, the Wigner function is positive for 
all times. Moreover, it peaks on the classical trajectories 
in phase space as the wave function spreads. The Wigner 
function can then be interpreted as a classical probability 
distribution for coordinates and momenta, showing sharp 
classical correlations at long times. In other words, the 
initial Gaussian state becomes highly squeezed and indis- 
tinguishable from a classical stochastic process. In this 
sense, one recovers a classical evolution of the inflaton 
rolling down the hill. 

A similar approach has been used by many authors to 
describe the appearance of classical inhomogeneities from 
quantum fluctuations in the inflationary era • Indeed, 
the Fourier modes of a massless free field in an expanding 
universe satisfy the linear equation 

0fc + (fc' - -)4>k - 0. (1) 

a 

For sufhciently long- wavelengths (fc^ <C a" /a), this equa- 
tion describes an unstable oscillator. If one considers an 
initial Gaussian wave function, it will remain Gaussian 
for all times, and it will spread with time. As with the 
toy model of Guth and Pi, one can show that classical 
correlations do appear, and that the Wigner function can 
again be interpreted as a classical probability distribution 
in phase space. (It is interesting to note that a similar 
mechanism can be invoked to explain the origin of a clas- 
sical, cosmological magnetic field from amplification of 
quantum fluctuations). 

However, classical correlations are only one aspect of 
classical behaviour. It was subsequently recognized that, 
in order to have a complete classical limit, the role of the 
environment is crucial, since its interaction with the sys - 
tem distinguishes the field basis as the pointer basis 
[We are reminded that, even for the fundamental prob- 
lem of the space-time metric becoming classical, simple 
arguments based on minisuperspace models suggest that 
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the classical treatment is only correct because of the in- 
teraction of the metric with other quantum degrees of 
freedom [l5l|.] 

While these linear instabilities cited above characterise 
free fields, the approach fails when interactions are taken 
into account. Indeed, as shown again in simple quantum 
mechanical models (e.g. the anharmonic inverted oscil- 
lator), an initially Gaussian wave function becomes non- 
Gaussian when evolved numerically with the Schrodinger 
equation. The Wigner function now develops negative 
parts, and its interpretation as a classical probability 
breaks down . One can always force the Gaussianity of 
the wave function by using a Gaussian variational wave 
function as an approximate solution of the Schrodinger 
equation, but this approximation deviates significantly 
from the exact solution as the wave function probes the 
non-linearities of the potential 0, 0| . 

When interactions are taken into account, classical be- 
haviour is recovered only for " open systems" , in which 
the unobservable degrees of freedom interact with their 
environment. When this interaction produces both a di- 
agonalization of the reduced density matrix and a posi- 
tive Wigner function, the quantum to classical transition 
is completed |l7| . 

In Ref.|^ we have considered an anharmonic inverted 
oscillator coupled to a high temperature environment. 
We showed that it becomes classical very quickly, even 
before the wave function probes the non-linearities of the 
potential. Being an early time event, the quantum to 
classical transition can now be studied perturbatively. In 
general, recoherence effects are not expected Tak- 
ing these facts into account, we have extended the ap- 
proach to field theory models 0, 13 • In field theory, one 
is usually interested in the long-wavelengths of the or- 
der parameter. Even the early universe is replete with 
fields of all sorts which comprise a rich environment, in 
the inflationary example, we considered a model in which 
the system-field interacts with the environment-field, in- 
cluding only its own short- wavelengths. This is enough 
during inflation. Assuming weak self-coupling constant 
(inflaton potential is flat) we have shown that decoher- 
ence is a short time event, shorter than the time tend, 
which is essentially the time by which diff'erent modes in 
the system sector cross the horizon. As a result, pertur- 
bative calculations are justified^. Subsequent dynamics 
can be described by a stochastic Langevin equation, the 
details of which are only known for early times |l9j . 

In our approach, the quantum to classical transition 
is defined by the diagonalization of the reduced density 
matrix. In phase transitions the separation between long 
and short-wavelengths is determined by their stability, 
which depends on the parameters of the potential. Dur- 
ing Inflation, this separation is set by the existence of the 
Hubble radius. Modes cross the aparent horizon during 
their evolution, and they are usually treated as classical. 
The main mativation of this talk is to present a formal 
way to understand this statement within the open qua- 
tum sytem approach. In the last sense, decoherence is the 



critical ingredient if we are to dynamically demonstrate 
the quantum-to-classical transition of the open system. 

The paper is organized as follows. In Section 2 we in- 
troduce our model. This is a theory containing a real 
system field (j), massless and minimally coupled to de Sit- 
ter background. We compute the influence functional 
by integrating out the environmental sector of the field, 
composed by the short-wavelength modes. Section 3 is 
dedicated to reviewing the evaluation of the master equa- 
tion and the diffusion coefficients which are relevant in 
order to study decoherence. In Section 4 we evaluate 
upper bounds on the decoherence times. As we will see, 
decoherence takes place before the end of the inflationary 
period. Section 5 contains our flnal remarks. 



II. THE INFLUENCE FUNCTIONAL AND THE 
DENSITY MATRIX 

Let us consider a massless quantum scalar field, mini- 
mally coupled to a de Sitter spacetime ds^ = a{r])[d'r]'^ — 
dx'^] (where ij is the conformal time dij — dt/a{t)), with 
a quartic self-interaction. The classical action is given by 



= / d^x a\rj) 



(2) 



where 0(77) = —l/lHrj) and (/)' — d(f>/dr] {a{rji) = \ [rji = 
— iJ^^], and H is the Hubble radius). Let us make a 
system-environment field splitting 



= 0<+0>, (3) 

where the system field contains the modes with wave- 
lengths longer than the critical value A~^ = 27r/Ac, while 
the bath field contains wavelengths shorter than A^^. As 
we set a{rji) = 1, a physical length Aphys = a('7)A coin- 
cides with the corresponding comoving length at the ini- 
tial time. Therefore, the splitting between system and 
environment gives a system sector constituted by all the 
modes with physical wavelengths shorter than the critical 
length Ac at the initial time rji. 

After splitting, the total action ||2Jl can be written as 



=5o[0<]+5o[</.>]+5int[</'<,0. 



(4) 



where So denotes the free field action and the interaction 
term is given by 

^i„t[0<,0>] = -A Jd^xaHv) {(l>%{x) + cl>i{x) 
+ 6<j>% {x)4>l {x) + A4>\ (a;)0> {x) + {x)cf>l (x) } (5) 

The total density matrix (for the system and bath 
fields) is defined by 

p[0+,0+,0<,(/);,t] = (0+0+|/5|0<0;), (6) 
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where and are the eigenstates of the field op- 
erators (/)< and 0>, respectively. For simplicity, we will 
assume that the interaction is turned on at the initial 
time rji and that, at this time, the system and the envi- 
ronment are not correlated (we ignore, for the moment, 
the physical consequences of such a choice, it has been 
disussed in |^). Therefore, the total density matrix can 
be written as the product of the density matrix operator 
for the system and for the bath 



(7) 



We will further assume that the initial state of the envi- 
ronment is the vacuum. 

We are interested in the influence of the environment 
on the evolution of the system. Therefore the reduced 
density matrix is the object of relevance. It is defined by 



Pr[0<><,'7] = y 250>p[0<,0>><>>,'7]- (8) 
The reduced density matrix evolves in time by means of 

d(f><i JrbhVi] pAm], (9) 



pAv] = J ' 



where Jr[?7,7yi] is the reduced evolution operator 



</ r'P<f 
i J4'<i 



X expl{5[0+]-5[0<]}F[0+,0-].(lO) 



The influence functional (or Feynman- Vernon functional) 



F[(f>~^, cj)^] is defined as 



5>/ 



>>/ r'p>f 



Vcp- exp(*{5[0+]+5i„t[0^,</'J]}) 

X exp (-i{S'[(/)>] -f 5i„t [(/)<, (?;'>]}). 

This functional takes into account the effect of the en- 
vironment on the system. The influence functional de- 
scribes the averaged effect of the environmental degrees 
of freedom on the system degrees of freedom to which 
they are coupled. With this functional, one can indetify 
a noise and dissipation kernel related by some kind of 
fluctuation-dissipation relation. This relation is impor- 
tant when one is interested in possible stationary states 
where a balance in eventually reached. During inflation 
we have a very flat potential well away from its minimum, 
and we are, in general, only interested in the dynamics 
over some relatively small time. For example, we would 
neglect dissipation during the slow roll period; but it is 
not correct during the eventual reheating phase. 

We define the influence action 0^] and the 

coarse grained effective action (CGEA) as 



F\ 



exp^M[(?!):t,0<], 



(11) 



(12) 



We will calculate the influence action perturbatively in A 
and we will consider only terms up to order and one 
loop in the h expansion. The influence action has the 
following form 

M[</.+ ,0-] = {(5i„t[0+ </>+]>o - (5i„t[0<,0;])o} 

+ ^{(^f„t[0^,0^]>o-[(5i„t[<^+,</.+])o]'} 

{{SintWt:, (/)>] 5*1111 [0<, 0>])o 

- (5int[0i,0^])o(5i„t[0<,0>])o} (13) 

i 



{'S'hit[0<>'?^>]>O - [{Sn,t[(t><,(t>>\)o\^^ 



where ( )o is the quantum average, assuming the envi- 
ronment is initially in its vacuum state. 

The influence functional can be computed, and the re- 
sult is 



ReM = j d'^^ (^^iv) P{x) 

+A2 j d^x J d^y a^{-q) a^{ri') 0{t] ~ 7?') 

X {64A3(a;)ReG^+(a;-2/)E3(y) 
+288A2(a:)ImGf+(a; - y)^2{y)} , (14) 

ImM = -A^ J d^x j d^y a^{rj) a\rj') 

{-32A3(a;)ImG^+(x-y)A3(2/) 

+ 144A2(x)ReGf+(x - y)A2{y)} , (15) 

where G^_^_{x — y) is the Feynmann propagator of the 
environment field, where the integration over momenta 
is restricted by the presence of the infrared cutoff A. We 
have also defined, 



p = 


^< 




A2 = 


^< 




S2 = 


^< 





1 



) ; A3 = ^ 
< ) ; S3 



<^) 



III. MASTER EQUATION AND DIFFUSION 
COEFFICIENTS 



In this Section we will obtain the evolution equation 
for the reduced density matrix (master equation), pay- 
ing particular attention to the diffusion terms, which 
are responsible for decoherence. We will closely follow 
the quantum Brownian motion (QBM) example |21ll22|. 
translated into quantum field theory 01 13- 
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The first step in the evaluation of the master equa- 
tion is the calculation of the density matrix propagator 
Jr from Eq. . In order to solve the functional integra- 
tion which defines the reduced propagator, we perform a 
saddle point approximation 

-'r[0</,0</,'7|0<i,0<i,?/»] ~ expiyl[(/)+^j,0<J, (16) 

where <p<ic\ solution of the equation of mo- 

tion SReA/S(l)'^\^+ — with boundary conditions 

^KciiVi) — (t>%i and </><ci('7) — '^</- Since we are working 
up to order, we can evaluate the influence functional 
using the solutions of the free field equations. This classi- 
cal equation is 0^ +27^0^ —V'^(/)< = 0, = a' (rj) / a{r])) . 
A Fourier mode 0^(a;) — J.j:. j^cp'i exp{ik.x}, satisfies 



a{rj) 



= 0, 



(17) 



where we have used V'^ = the fact that 

^^j^^ — It is important to note that for longwave- 
length modes, fc ^ 2/77^, Eq. lfT7|) describes an unstable 
(upside-down) harmonic oscillator (T3 |. 
As a classical solution can be written as 



•^f W = 0f(fc)ui (77,77/) +4(fc)u2(??,r;y), (18) 



These identities allow us to remove the initial field con- 
figurations (j>f' , by expressing them in terms of the final 
fields (j)f and the derivatives , and obtain the master 
equation. 

The full equation is very complicated and, as for 
quantum Brownian motion, it depends on the system- 
environment coupling. In what follows we will compute 
the diffusion coefficients for the different couplings de- 
scribed in the previous section. As we are solely inter- 
ested in decoherence, it is sufficient to calculate the cor- 
rection to the usual unitary evolution coming from the 
noise kernels (imaginary part of the infiuence action). 
The result reads 

ihdnPM<f^'l><f^'n\^ {<l><f\ Hrcn,Pr |0<f) 



Fi Di(fco;??) + r2 D2{ka\ii) 



(21) 



where we have defined Fi = 



k-3 \2 



IP' 



and r2 = 



where 



. V is the spatial volume inside which 

there are no coherent superpositions of macroscopically 
distinguishable states for the system field. The ellipsis 
denotes other terms coming from the time derivative that 
not contribute to the diffusive effects. This equation con- 
tains time-dependent diffusion coefficients Di{t). Up to 
one loop, only Di and D2 survive. Coefficient Di is re- 
lated to the interaction term (/><i^>, while D2 to 
These coefficients can be (formally) written as 



sm[k{T]~r]f)](l + krir]f)+coa[k{r]~rif)]{T]f-r]) Di{ko,T]) = 2-j / dr,'a'^{r])a*{Tj')F^i{ko,T],7]') 



Ml 



U2 



Idem num. 7/ — > Tyi 
sin[A:(77i - 77)] (^ + krjrji) + cos[k{r] - 77^)] (77 - 7;^) 



Idem num. 77 



We will assume that the system-field contains only one 
Fourier mode with k = kg. This is a sort of "minisu- 
perspace" approximation for the system-field that will 
greatly simplify the calculations, therefore we assume 



and 



X ImG^^+(3fco,7;,7?') ^(3fco-A), (22) 



D2iko,rj) = 36^ J\v'a\ri)a\rj')F^iiko,v,ri') 

X [RcGf+(2fco,7?,^') + 2RcGf+(0, 77,77')] , (23) 
where the function Fc\ is 



(/.±^^(f,77) = 0f' (77)cos(^o.x), (19) 

where (j)^^^ is given by H18|) . 

In order to obtain the master equation we must com- 
pute the final time derivative of the propagator Jr. After 
that, all the dependence on the initial field configura- 
tions (coming from the classical solutions 0<'^') must 
be eliminated. In previous publicatiuons, we have shown 
that the free propagator satisfies Q 



Fci{ko,r],r] ) ^ h — cos /co(77 - 7/ ) . (24) 

koV V 

The explicit expression of these coefficient are compli- 
cated functions of conformal time, the particular mode 
fco, and the cuto ff A, and we will show them in a sep- 
arate publication z3j. For the porpuse of this talk, we 
will use some analytical approximations, which will allow 
to obtain an estimation of the scale of decoherence in a 
particular case. 



m' 1(77/, 77/) 
Jo- 



"i('7,^/)] 

(20) 



IV. DECOHERENCE 

The effect of the diffusion coefficient on the decoher- 
ence process can be seen considering the following ap- 



5 



proximate solution to the master equation 



X exp 



dr/ Dj{kQ,A,r]) 



(25) 



where p" is the solution of the unitary part of the master 
equation (i.e. without environment), and Tj includes the 
coefficients in front each diffusion term in Eq.|2lJ. The 
system will decohere when the non-diagonal elements of 
the reduced density matrix are much smaller than the 
diagonal ones. 

The decoherence time-scale sets the time after which 
we have a classical field configuration, and it can be de- 
fined as the solution to 



'to 



dr]Dj{ko,A,r]). 



(26) 



these equations are obtained under the following condi- 
tions 



%i U' 



2 U" , 



(30) 



where = 1/G is the Plank mass. 

Definig the end of the inflationary period setting ejj ^ 
f , one can set (f>{N) w VSTV mpi, where N = \na[rif)/ a{rf) 
is the e-fold number. Thus, we assume the mean value of 
the system field at time of decoherence is </>, and we set 
((/)+ - (^<) ~ 10-^(/), and V ~ H-^. 

From previous considerations, we can show 



^""^ - 6i/ 1 A21O-1506 



and 



6H 



In- 
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, (31) 



We will solve this equation only in a very approximated 
way in order to find upper bounds for the decoherence 
times coming from each diffusion coefficients. A more 
refined evaluation of decoherence times will be shown in 
Ref. |2^. For example, when we consider a mode fco < 
-ff, we can probe that a good approximation to Di is 
given by 



i^r (fco~i?,^) 



(1 + Htj) 



(27) 



This is valid for modes shorter or the same order than H] 
but it is an overestimation for k^ <C H . As Di is defined 
under the constraint A/3 < k^ < K, and considering we 
will set the critical length A < H, ours is a good estima- 
tion. Analysing the coefficient D2 , it is possible to find an 
approximated expression for small values of fco (respect 
to H). For ko/H > 1, that corresponds to those modes 
out of the aparent horizon at rji, the diffusion coefficient 
D2 is an oscilatory function and it has a maximun when 
fco ~ A, this is the behavior noted for conformally cou- 
pled fields |2^. Finally, for long- wavelength modes, we 
can write 



tD2 < — In 



3H 




(32) 



Using N — Htcnd > 60 as an estimative scale to the 
end of inflationary period. 



t, 



end 



< ln{AfOi^} < f , 



(33) 



corresponding to the diffusion term Di and values of A 
10^^. The scale coming from D2 is 



tD2 ^ J_ 
tend ~ 10' 



(34) 



From scales , and tn^ , we can see that for a given 
mode ko < A < H, decoherence is effective by the time 
in which inflation is ending. 



FINAL REMARKS 



Dr-iko<H,,)^j^. (28) 

Both approximations above, are close to the exact co- 
efficients when A < H. 

In order to quantify decoherence time we have to fix 
the values of Fi 2- For this, we have to assume values to 
A, V, A, and E. We will use the more conservative choice 
in order to have a lower bound to the decoherence time. 

Assuming slow roll condition l/2{d<j)/dt) <^ U = A(/)^, 
the classical equations (using cj) <^ U') are 



Let us summarize the results contained in this paper. 
After the integration of the high frequency modes in Sec- 
tion 2, we obtained the CGEA for the low energy modes. 
From the imaginary part of the CGEA we obtained, in 
Section 3, the diffusion coefficients of the master equa- 
tion. System and environment are two sectors of a single 
scalar field, and the results depend on the "size" of these 
sectors, which is fixed by the critical wavelength A~^. 

In Section 4 we analysed the decoherence times for 
those modes in the system whose wavelength is shorter 
than the critical value. 

We have shown some analitycal approximations that 
allow to conclude that if we consider a critical length 
A H, those modes with wavelength fco <C A are the 
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more affected by diffusion throught coefficient D2- For 
these modes, we can show that the effect is not dependent 
of the critical A pl |. 

If one consider a cutoff A > H, and modes H < ko < 
A, diffusive effects are larger for those modes in the sys- 
tem whose wavelength is close to the critical A~^ [23 |. 

In Ref. we will present the complete expression 
of the diffusive terms, an also an extensive analysis of 
the evaluation of the timescale for decoherence in several 
cases of interest. 
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